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not understand crystals very deeply. 
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In the last lecture, I explained how condensed-matter 
and high-energy physicists used topological theories to 
describe defects excitations in solids. In this lecture, I'm 
going to make fun of topology. Q Actually, I'm going to 
start by talking about constraints, then "massive" fields 
and how they produce constraints. I'll then turn to the 
Meissner-Higgs effect in superconductors, and finally ex- 
plain why I don't understand crystals. 



I. CONSTRAINTS 

Consider figure 1. See the beautiful ellipses and hyper- 
bolas? Remember that topology treats ellipses as rubber 
bands. Any topological theory has got to miss the key 
feature of the beautiful structures produced here: the ge- 
ometrically perfect ellipses with dark lines coming out of 
one focus. 

Figure 1 is a photograph of a drop of fluid, squeezed 
between two microscope slides. The microscope is fo- 
cused, let's say, on the surface between the fluid and the 
bottom microscope slide: the ellipses are stuck onto the 
glass. The sizes of the ellipses are roughly given by the 
thickness of the fluid layer. The fluid is a smectic A liq- 
uid crystal. deGennesQ has a fine discussion and some 
nice pictures too. 

In 1910, Friedel figured out why this liquid crystal 
forms these geometrical structures. He learned all he 
needed to know from his high-school geometry class. He 
actually worked backward, and used the ellipses to de- 
duce what kind of broken symmetry the liquid had. Since 
none of you were taught about the cyclides of Dupin in 
high school[3, I'd better start with the broken symmetry 
and work forward. 

Smectic liquids form equally spaced layers. Some of 




FIG. 1: Ellipses: Defects in a Liquid Crystal. This is a drop 
of smectic A liquid crystal, squeezed between two microscope 
slides. The microscope is focused on the surface of the drop, 
where it contacts the glass. Notice the beautiful, geometrical 
ellipses. Notice that a line seems to exit from the focus of each 
ellipse. This line turns out to be a hyperbola (figure 4). The 
visible ellipses and the hyperbolas are where the smectic layers 
pinch off to form cusps. These defects are not topological: 
they are geometrical consequences of the constraint of equal 
layer spacing. From |^, figure 7.2, photo by C. Williams. 



them are compounds that, like soap, naturally form mem- 
branes and films: I think smectic is the Greek word for 
soap. Others are long thin molecules like nematics, which 
for some reason not only line up, but segregate into planes 
(figure 2). The molecules have liquid-like order in the 
planes. Like crystals, they have a broken translational 
symmetry, but only in one of the three directions. 

Now, the important excitations for smectics are those 
that bend the layers. In figure 3, we see a two- 
dimensional analogue of the smectic liquid crystals: 
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FIG. 2: Order in Smectic Liquid Crystals. Smectic liq- 
uid crystals are formed of layers of molecules. In each layer, 
the molecules are in a random, liquid configuration. Crystals 
have broken translational symmetry along three independent 
axes: smectic A liquid crystals have broken translational sym- 
metry in only one direction (normal to the layers). 

equally spaced curves in the plane. Suppose we start 
with one curve and work outward. As you can see from 
the figure, the next curve is not precisely the same shape: 
keeping the surfaces at an equal spacing makes concave 
regions become sharper and convex regions become more 
rounded. It is easy to see that eventually the concave re- 
gions will become pinched: these pinches are the defects. 
They are not topological defects, since rounding them a 
bit makes them go away: they are geometrical defects 
produced by the constraint of equal layer spacing. 

Most curves, like the one shown in figure 3, form one- 
dimensional pinched regions: only concentric circles and 
structures made from them can keep the pinched regions 
to points. In three dimensions, the only equally spaced 
surfaces with points as pinched regions are concentric 
spheres. Now, what Friedel knew and you don't know 
is that the only 3-D surfaces with one-dimensional line- 
like defects are the cyclides of Dupin, |fl and the pinched 
regions form ellipses and hyperbolas. |lf| 

Figure 4 shows the cyclides of Dupin. Notice that they 
pinch off on two curves: an ellipse and a hyperbola. The 
hyperbola is perpendicular to the plane of the ellipse, and 
passes through its focus. That's what you see streaming 
out of the foci in the photo, and why you don't see one for 
each focus. My contribution to the field (with Maurice 
Kleman) was to realize that these cyclides of Dupin fit to- 
gether nicely inside concentric spheres, which explained 
neatly the ways the ellipses always seemed to fit together 
(figure 5). Maybe the concentric spheres form because 
the layers nucleate on a dust particle on one of the mi- 
croscope slides: when the spheres touch the other slide, 
the concentric spheres get twisted (they like to sit per- 
pendicular to the glass) and the ellipses and hyperbolas 
form to relieve the strain. 

Now, why do I show you this? It isn't just to show that 
there is more to the world than topology. Mostly, it's to 
illustrate the two themes of this lecture: constraints, and 
expulsion. 

If we define an order parameter n for the smectic to be 
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FIG. 3; Equally Spaced Layers: Defect Formation. 

Here we consider a two-dimensional analogue of a smectic 
liquid crystal. The smectic layers are represented by curves 
in the plane. The lowest energy state, of course, consists of 
parallel straight layers, but the layers often settle into more 
complicated patterns, with defects. For reasons that we dis- 
cuss in this lecture, and which are not completely understood, 
smectic layers will deform by bending, but will remain strictly 
equally spaced (except very near boundaries and defects). 
The constraint of equal layer spacing has weird nonlocal con- 
sequences. First, one can see that as one moves outward the 
concave regions become more pinched, and eventually form 
cusps. Second, one can see that a line perpendicular to one 
layer (a generator) will be perpendicular to the next one too. 
These generators intersect on a surface known as the evolute, 
and it is when the layers hit the evolute that a defect oc- 
curs. As one sees here, the defect is a line of pinched surfaces: 
in three dimensions it is typically a two-dimensional surface. 
This costs lots of energy. The only way in two dimensions 
to have a point-like low-energy defect is to have concentric 
circles: only circles have zero-dimensional evolutes. The only 
way in three dimensions to have one-dimensional evolutes |3] 
is to have cyclides of Dupin: the defects are ellipses and hy- 
perbolas passing through one another's foci (figures 1 and 4). 

the unit normal to the smectic layers (n^ = l)j then the 
constraint that the layers be equally spaced implies 

(duz/dy — duy/dz \ 
dux/dz — duz/dx =0. (1) 
duy/dx - dn^/dy J 

(This is derived, for those who know a bit about vector 
calculus, in the appendix.) This is a remarkably pow- 
erful constraint. For example, knowing the position of 
one layer determines all the others! We show this math- 
ematically in the appendix, but you saw it physically in 
figure 3: given one layer, there is only one way to place 
the next one preserving exactly equal spacing. 

There is a pretty good analogy here to analytic con- 
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FIG. 4: Focal Conic Defect. Here we see the smectic sur- 
faces which form the focal conic defects seen in figure 1. These 
are the cyclides of Dupin. The surfaces go from banana- 
shaped to squashed doughnuts to apple-shaped. The points 
on the bananas and the dimples at the stem and bottom of 
the apples are defects, which scatter light and show up in 
figure 1. (Only the dimples of the apple are shown.) The 
banana defects lie on an ellipse, and the apple defects lie on 
a hyperbola which passes through the focus of the ellipse. 
Usually, the whole pattern isn't found in the experimental 
sample. As you see in figure 1, the domains aggregate to- 
gether in clumps. Each ellipse in figure 1 has a conical region 
for its smectic layers. 

tinuation. For those of you who know about complex 
analysis, you know that an analytic function obeys the 
Cauchy-Riemann equations. If we let n{x + iy) = 
nx{x + iy) + iny{x + iy) be an analytic function, then 

/ dn^/dx- dny/dy\ 
\dnx/dy + duy/dx ) 

As you know, analytic functions have really bizarre prop- 
erties. If you know an analytic function in a small region, 
you can figure it out everywhere else, just like the order 
parameter in smectics. The point singularities of ana- 
lytic functions have a rich and interesting classification 
(simple poles, essential singularities, ...). Both in ana- 
lytic functions and in our smectic problem, constraints 
on the derivatives of our order parameters produced re- 
ally bizarre, nonlocal, geometrical consequences. 

II. MASSIVE FIELDS 

We've discovered that constraints can have beautiful, 
geometrical consequences. How are the constraints en- 
forced? Clearly, it is possible to stretch the smectic lay- 
ers apart, or to compress them together: why doesn't 
this happen in practise, especially when the layers are 
being bent and twisted? The curl of h is constrained to 
zero. Why are magnetic fields pushed completely out of 
superconductors? The magnetic field is constrained to 




FIG. 5: Focal Conic Defect Meshing Onto Concen- 
tric Spheres. The conical regions in figure 4 combine into 
compound defects by meshing onto the concentric sphere de- 
fect. Concentric spheres axe the only surfaces with zero- 
dimensional defects. The surfaces on the edges of the cones 
mesh smoothly onto the concentric spheres. 

zero. Why isn't it possible to find an isolated quark in 
nature? Quarks have non-zero "color" , and the net color 
is constrained to zero. 

These constraints come from minimizing the energy. 
Saying that magnetic fields can happen inside supercon- 
ductors is just like saying that marbles can sit on the 
side of a hill: it can happen, but not if the marbles are 
allowed to roll to minimize their energy. Under what 
conditions does the energy enforce a constraint? We say 
that it happens when the order parameter field develops 
a wMss. We'll explain this term in a moment, but let's 
first give a simple example. 

Suppose we have a fluid in one dimension. The density 
of a fluid is the important variable in describing its state. 
Suppose the density of the fluid is po + p{x), where po 
is the ideal density (which the fluid would have if left to 
itself) and the order parameter p{x) describes the devia- 
tion from the ideal density. A sensible free energy might 
be 

Sfiuid = Jdx (l/2)(dp/dx)2 + (l/2)mp2. (3) 

The flrst term in the energy resists sudden changes in the 
density: having a high density region right next to a low 
density region costs extra. The second term in the energy 
says that deviations from the mean density cost energy, 
with m a coefficient which says how much deviations cost. 
Unlike phonons, where the order parameter u{x) could 
be uniformly shifted without energy cost, here the lowest 
energy state happens when the density is at its mean 
value p{x) = 0. 
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What happens when we try to find the minimum en- 
ergy state? Clearly the best we can get is the ideal state 
p{x) = 0, which has zero energy £ fluid- Perhaps, though, 
we're pulling on the density at the two ends (figure 6). 
If the liquid is in a trough of length L, we'll insist that 
p{Q) = Pi and p[L) — pj. What configuration p{x) mini- 
mizes the energy then? Clearly, it should sag towards pa 
inside, but how? 




FIG. 6: Massive Fields Decay Exponentially. Minimiz- 
ing the energy Ejiuid in equation|21 with boundary conditions 
p(0) = pi and p{L) = Pf- It is easy to understand physically 
what is happening. The system wants to achieve p = 0, and 
it sags to that value as quickly as it can, balancing the costs 
of (dp/dx)^ energy against the gain. The solution decays ex- 
ponentially to zero with a decay constant y^m. 



Here I'll show you a simple case of what's called the 
calculus of variations. I apologize for the math, but it 
is really a useful method. The trick is to realize that if 
p{x) is the minimum energy configuration, then p(x) + 
d{x) must have a higher energy, whatever S{x) we might 
choose. 



from the interior, pulling it on the boundary only affects 
a region of length y/m, and the order parameter expo- 
nentially decays into the bulk, p is constrained to zero 
in the inside of the sample! 

Why do we call this a mass? The name comes from 
particle physics. The photon is massless. Two charges 
ei and 62 separated by a distance r interact by a force 
whose magnitude goes as eie2/r^: this is Coulomb's law. 
The particle physicists interpret this force in terms of the 
two particles exchanging "virtual" photons. (I think of 
the decay as the virtual photons being diluted over a 
sphere of radius r.) Now, the strong interaction between 
protons an neutrons has a different form: the force be- 
tween them is always attractive, and goes as e~^''/r^. 
The exponential decay is extremely important, since it 
keeps the nuclei of different atoms from attracting one 
another. (We'd all have collapsed into neutron stars or 
worse were it not there!) At long distances, the particle 
physicists interpret this force as the proton and neutron 
exchanging virtual Dions.[lll| Since the pion isn't mass- 
less, the virtual pion field decays exponentially for ex- 
actly the same reason that p(x) decayed in our example 
above. 

So, to enforce a constraint, we need to give the corre- 
sponding field a mass. Let's see how that is done. 



III. THE MEISSNER-HIGGS EFFECT 



Sip + 6) - Sip) (4) 
— I idp/ dxd5 / dx -\- mpix)5{x) 



+ il/2)id5ldxf + (l/2)m(52) dx 
> 0. 

Now, if we confine our attention to small Six), we can 
ignore the last two terms (because they are quadratic, 
rather than linear, in 5). The first term we integrate by 
parts, so 



L L 

dx dS/dx dp/dx — iSdp/dx) 



dx S d^p/dx^ 



(5) 

Now, S mustn't change the values at the endpoints, so 
(5(0) = (5(L) = and the boundary terms inOdrop out. 
We're left, then, with the equation 

Sip+S)-S{p) [ dxi-d^p/dx^+mpix))Six) >0. (6) 



Now, this must be true for any Six) we choose. This can 
only happen if —d^pjdx"^ + mpix) ~ 0, so p" ~ mp. 

The solutions to this equation are, of course, p — 
j^^-^/mx _|_ j^gVmx^ -y^g Vary the arbitrary constants 

A and B to match the boundary conditions p(0) = pi 
and piL) = Pf, and we see (figure 6) that p is expelled 
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FIG. 7: Superconductors Expel Magnetic Fields. A 

magnetic field passing through a metal will be pushed out 
when the metal is cooled through its superconducting tran- 
sitions temperature. This can happen in two difi'erent ways. 
In type I superconductors like lead (chemical symbol Pb), the 
superconductivity is pushed entirely outside the sample. In 
type II superconductors like niobium (Nb), the magnetic field 
is broken up and confined to defect lines called vortices. In 
both cases, the magnetic field is swept out of the remain- 
der of the sample. The magnetic field penetrates a distance 
A ~ lOOAinto the sample from the boundaries or from the 
vortex lines. 

In this section, I want to explain how superconduc- 
tors expel magnetic field. This is a really beautiful 
argument, which I've basically taken from Coleman's 
presentation. I'm afraid that there is some math and a 
lot of physics that I need to introduce. Most of you will 
get lost, but the pictures will be nice anyhow. 
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A. Introduction to the Meissner Effect. 

Superconductors are named for their ability to carry 
currents of electricity with absolutely no losses. They 
have another, closely related property which is no less 
amazing: they are a perfect shield for magnetic fields. 
Remember the old science fiction stories about the sci- 
entist who finds a material which is impervious to the 
gravitational field, paints the bottom of his spacecraft 
with it, and falls to the moon? Superconductors work 
that way for magnetic fields. 

Ashcroft and Mermin have a nice, not too technical dis- 
cussion of superconductors in one of the last chapters in 
their textbook. Q Figure 7 shows the two types of super- 
conductors, represented by lead and niobium. At high 
temperatures, when the materials aren't superconduct- 
ing, the magnetic field penetrates the materials almost 
as if they weren't there. (Iron would pull the magnetic 
field lines inward.) Lead, when superconducting, pushes 
the magnetic field out: just as for the example in sec- 
tion II, the field a distance r inward from the boundary 
decays like B = BQe~^^^. If you put too high a field, 
the lead will give up and let the field in: but it will stop 
superconducting. 

On the right, we see that niobium behaves a bit differ- 
ently. It expels small magnetic fields like lead does, but 
larger fields are pushed into thin threads, called vortex 
lines. These two general categories are (rather unimagi- 
natively) called type I and type II superconductors. The 
vortex lines are the topological defects for the supercon- 
ductor (lecture 1). Superconductors are described by a 
complex number tp — pe'^, whose magnitude — p is 
roughly constant. The order parameter at low tempera- 
tures is the phase 9, and thus the order parameter space 
is a circle . A vortex line must pass through any loop 
around which the phase of the order parameter changes 
by 27r. The magnetic field in type II superconductors de- 
cays like B = Boe~^^^ where here r is the distance to the 
vortex line. Magnetic field is squeezed out of the bulk of 
the material into these defects. 

So, the magnetic field isn't actually stopped, it just 
peters out. What kind of a leaky shield is that? Actu- 
ally, it's about as good as one can hope: after all, the 
magnetic field won't be able to tell it's in a supercon- 
ductor until it gets inside a bit! (Atoms don't go super- 
conducting, only huge heaps of atoms together, so the 
field has to pass through a heap or two to realize that 
it isn't wanted.) Anyhow, A is usually pretty small, a 
few hundred Angstroms or so. An 0.1mm thin layer of 
superconducting paint naively would let through a field 
one part in e-'^°°°° ~ 10-"°*^° of the original. Unfortu- 
nately, it usually doesn't work so well: a few vortex lines 
get stuck on junk in the paint, and let in comparatively 
large fields. 

Before we can explain the repulsion of magnetic fields, 
we should explore the broken symmetry. Let's start with 
superfluids, which are simpler. 



B. Superfluid Free Energy and Spontaneous Symmetry 
Breaking. 




B 

FIG. 8: Superfluid Free Energy, (a) T > Tc. Unbro- 
ken Symmetry. The free energy for a normal metal or fluid, 
above the superconducting or superfluid transition tempera- 
ture, for a uniform order parameter field -i/;. The vertical axis 
represents the energy alipl^ -\- and the horizontal axes 

represent the real and imaginary parts of 1/). The coefficient 
a > 0, so the minimum of the energy is at = 0. Notice that 
the energy is invariant under the symmetry ■;/> — > e^^tp (cor- 
responding to rotating the figure about the vertical). This is 
a symmetry of the free energy. Notice also that the lowest 
energy state t/; = is also unchanged by this rotation: the 
symmetry is unbroken above Tc- 

(b) T < Tc- Broken Symmetry. The free energy E super fluid 
for helium below the superfluid transition temperature. The 
energy now looks like a Mexican hat: it is still invariant under 
rotations about the vertical axis. Since now a < 0, the energy 
is at a minimum along a circle, of radius \ip\ — a/2/3 and 
arbitrary phase 6. The superfluid must choose between these 
various possible phases, and that choice breaks the symmetry. 
This is a good example of spontaneous symmetry breaking: 
just as the magnetization of a magnet selects a direction in 
space and breaks rotational invariance, the superconductor 
picks out a value of 9. 

The order parameter for a superfluid, just as for a su- 
perconductor, is a complex number ip. The free energy 
for the superfluid is usually written as|l2| 

Ssuperflmd = / dV\V^'' + a|V'|' + ■ (7) 
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Above the superconducting transition temperature Tc, 
the coefficient a > 0. If we imagine a constant order 
parameter field, the free energy forms a bowl (figure 8a) 
with a minimum at zero, as a function of the real and 
imaginary part of ip. Zero order parameter corresponds 
to a normal metal (for a superconductor), or a normal 
liquid (for a superfluid). 

Below Tc, a < 0, and the potential is at a minimum 
for po — \tp\ — ■y/a/2/3: the potential in the complex 
plane looks like a Mexican hat (figure 8b). Now there 
are many possible ground states: for any 9, a constant 
order parameter field ijj = Poe'^ is a ground state. Be- 
cause the free energy depends only on [ipl and |V';/'|, it 
is symmetric to changing the phase 9: the superconduct- 
ing state chooses a specific value for 9, and thus sponta- 
neously breaks the symmetry. The circle of ground states 
in the brim of the Mexican hat is the order parameter 
space for the superconductor. 

We can write the free energy in terms of 9: 

Ssuperflmd = J dV\V + p^\V9\^ + ap" + Pp\ (8) 

As we discussed in the previous section, p is "massive" . 
In figure 8b, if we vary p slightly away from po, the en- 
ergy increases quadratically: ap^ + fip"^ — {apQ + (3pq) « 
(a -I- &(5pq){p — Po)^ The effective free energy for p near 
Po is precisely of the formal (except for unimportant con- 
stant shifts), with TO = a -f 6/3pQ. Thus just as before, 
p will rapidly be drawn to its minimum energy state po- 
Because p is massive, it is basically constrained to stay 
at its minimum value. This is why it is ignored at low 
temperatures in writing the order parameter field. Here, 
the constraint doesn't do anything interesting: our next 
constraint will be more interesting. 

The 9 field keeps the symmetry of the original model: 
rotating \t to 9 + 9q doesn't change the energy a bit. 
It is a Goldstone mode for our problem, and long- 
wavelength plane waves produce what is known as "sec- 
ond sound" in superfiuids. Second sound turns out to be 
heat waves: pulses of temperature which propogate like 
waves through the superfluid. 

C. Superconducting Free Energy and the Higgs Mecha- 
nism 

To describe the expulsion of magnetic field from super- 
conductors, I have to tell you how magnetic fields interact 
with the superconducting order. I'm afraid this will be 
rather sketchy, and I apologize for trying. 

First of all, the particles which superconduct are pairs 
of electrons. Electrons are charged, and repel one another 
with electric fields. Thus the electrons interact with elec- 
tric fields. We learn in the second semester of physics (if 
we're lucky) that electric and magnetic fields are closely 
related to one another. (This was discovered by Einstein: 
a moving electric E field develops a magnetic B compo- 
nent.) 



Now, the E and B fields can be written at the same 
time in terms of another field A. It is this new field which 
is easiest to work with. In particular, 

B = curl A (9) 
^ ^dA^_dAy dA^_dA^ dAy dA,^ 
dy dz ' 9z dx ^ dx dy 

The magnetic energy is Smagnetic = J dV B'^. 

Now, you remember that I mentioned earlier that light 
(photons) is massless? You may know that light is some- 
times called "electromagnetic radiation" . The "order pa- 
rameter field" for light is precisely the A field. We can 
see by expanding in terms of A that the energy for 
the A field 

doesn't have any terms like A^. When we add the energy 
from the electric fields, this is still true: light is massless 
because the electromagnetic energy involves only deriva- 
tives of A. 

Now, I need to know how the electromagnetic order 
parameter A interacts with the superconducting order 
parameter ip. I'll just tell you. The free energy for a 
superconductor looks like 

S superconductor = J dV \\7 - iA^P]"^ + a\i;\^ + Pl^P]* + B^ 

(11) 

If we set ip = 0, we get the magnetic energy B^ for the 
A-field. If we set A = 0, we get the superfiuid energy 
[7| I don't know of a way to motivate the way in which 
we couple the A field to the gradient VV'- I don't think 
anyone has a simple derivation. This way of connecting 
the two is called "minimal coupling" , which just gives a 
name to the unexplained fact that the simplest way of 
coupling the two gives the right answer. 

Now, if we assume T < so a < and p ^ Poe'^, we 
find 

j dV pl\V9-A\'^ + {dA,Jdy~dAy/dzf+---. (12) 

We want to know if A or 6* is going to develop a mass. 
The problem is, E superconductor doesn't look quite like the 
form 121 for either one. If we combine the two into a new 
order parameter field C — V9 — A, and use the fact that 
the second partial derivative d^9/dzdy = d^9/dydz, we 
see that 

curie = {dC,/dy-dCy/dz, • • •) 
= [dAJdy-dAy/dz, ■■■) 
= B (13) 

so 

f « [ dV pIC^ + {dC/dy - dCy/dzf + ■ ■ ■ . (14) 
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Thus the new, combined field C is massive. C will be 
constrained to zero in the bulk, exponentially decaying 
like CqCP"^. The magnetic field B = curlC thus also 
decays, and the penetration depth A = l/po- 

We started with a massless photon field A and a mass- 
less Goldstone mode 9. We ended up with only one field 
C, with a mass. Did we lose something? No, actually 
C has three components: two components corresponding 
to the original two polarizations of light, and one compo- 
nent corresponding to the Goldstone mode. Coleman*^ 
says "the Goldstone boson eats the photon, and gains a 
mass" ! 

The Weinberg-Salaam theory of the weak interaction 
is exactly analogous to the theory of superconductivity. 
The role of lead or niobium is played by the vacuum. 
The free energy of the universe has an SU{S) symmetry, 
which is spontaneously broken to a smaller symmetry 
SU{2) X C/(l). The W"^ and Z bosons which now mediate 
the weak interaction used to be massless: they and the 
photon were all part of one big ^-field. If current theories 
of cosmology are true, this "superconducting" transition 
occurred in the first instants after the Big Bang. 

Now, after explaining superconductors, the weak in- 
teraction, and the phase transition in the early universe, 
let's return to why we don't understand crystals. 

IV. THE MYSTERY OF THE CRYSTALS 




FIG. 9: Polycrystal. Many crystalline materials, such as 
metals, normally aren't made of a single crystal. They are 
formed from many crystalline domains: a poly crystalline con- 
figuration. I show a schematic of a polycrystal here. The 
important thing to notice is that the atoms within a domain 
are almost undeformed except right next to the domain wall. 
All the rotational deformation is expelled into sharp domain 
boundaries. 

Normally, when you think of crystals, you think of dia- 



monds, snowflakes, or maybe salt crystals. [I3 These are 
single crystals: the sodium and chlorine atoms in a grain 
of salt sit in registry all the way across the grain, giving 
it its cubical shape. Did you know that metals form crys- 
tals? In the last lecture, I mentioned dislocation lines in 
a copper crystal. Metals don't have big facets and cor- 
ners because they are polycrystalline. The atoms in a 
metal also sit on a regular lattice, but the metal breaks 
up into domains in which the lattices sit at various an- 
gles (figure 9). Because there are lots of small domains, 
copper doesn't form facets like salt grains and snowflakes 
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FIG. 10: Growing a crystal from a liquid: forming a 
polycrystal. Polycrystals can form for lots of reasons. If 
one cools a liquid quickly, one can find that crystalline re- 
gions can form in many different places almost simultaneously. 
Since they will have random orientations, they won't match 
up when they meet. When they do meet, rearrangements 
of atoms will occur to try to realign and merge the domains 
(coarsening). As we continue to cool and wait, this process 
will eventually stop, leaving us with different domains. 

What Ming Huang (one of my students ||7|) and I have 
been trying to explain for years is why those little do- 
mains form. It's easy to see that different regions might 
grow with different orientations (figure 10). When they 
touch, the different domains will start pushing and twist- 
ing one another, trying to make one big domain. It isn't 
hard to believe that they will stop growing after a while, 
fighting one another to a standstill. What we've been 
trying to understand, though, is why the final state is 
made of perfect little crystals separated by sharp domain 
walls. 

Now, I don't want to exaggerate. There are perfectly 
good explanations for why crystals form domain walls. 
They just aren't as beautiful and general as they might 
be. They don't fit in with the general ideas of broken 
symmetries and order parameters: they apply only to 
crystals. Our explanation for why superconducters don't 



have a Goldstone mode was perfectly OK before Higgs 
came too. He made it beautiful and generalized it to ex- 
plain something completely different. Ming and I want 
to understand grain boundaries in a way which will make 
simple and clear where else similar phenomena might oc- 
cur. At least, we'd like to understand why focal conies 
occur at the same time. Domains formed by breaking 
translational symmetry in one direction and in three di- 
rections should have the same kind of explanation! 




FIG. 11: Domain Wall. Here we see a single domain wall. 
Notice that the domain wall can be also thought of as a series 
of dislocations. The strain field inside the crystal due to a 
Une of dislocations can be shown to decay exponentially, just 
as the magnetic field dies away around a vortex line. 
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FIG. 12; (a) Rotational distortion of a crystal. If we 

take a thick piece of metal and rotate one end with respect to 
another, it will start by bending uniformly. As it continues to 
bend, dislocations will form to ease the bending strain. These 
line dislocations will start off distributed irregularly through 
the sample, (b) Domain walls form to expel rotations. 
If we hold the rotation for a long time, and let the dislocations 
move around, they will lower their energy by arranging them- 
selves into domain walls. Between the domain walls we find 
undistorted crystal. This process is called polygonalization. 



Figure 11 shows a domain wall in a crystal. The crys- 
talline ground state rotates as one crosses the domain 
wall. The atoms at the wall are quite unhappy. You'd 
think that they would push and pull on their neighbors, 
and that there would be strains leaking far into the crys- 
tal. This isn't true. In fact, there is a well-known rule in 
the materials science literature, that the strain field from 
a domain wall dies away exponentially as one enters the 
grain. 

Doesn't that sound like a Meissner effect? 

There are more analogies. Crystals break both the 
translational and the rotational symmetry of liquids. 
Many liquid crystals only break the rotational symme- 
try. They have Goldstone rotational waves: if you rotate 
a large region inside a liquid crystal, it will cost little en- 
ergy, and will slowly rotate back. When the translational 
symmetry is also broken, the rotational Goldstone mode 
disappears! If I rotate one piece of a crystal with respect 
to another, it costs an enormous energy (figure 12a). If I 
let the distorted crystal rearrange locally to reach equi- 
librium, the rotational deformation is expelled into grain 
boundaries (figure 12b), a process known in the field as 
polygonalization. Just like the massless photon devel- 
oped a mass when the superconducting transition broke 
the gauge symmetry, the massless rotational mode devel- 
ops a mass when the translational symmetry is broken. 



This is surely also related to some of the old problems 
in the topological theory of defects. In describing a crys- 
tal, everybody uses the displacement field u(x) and its 
derivatives. Now, as we saw in lecture 1, u{x) describes 
the broken translational order, but not the broken oricn- 
tational order. Why don't we also have a rotation ma- 
trix R(x)7 For example, in figure 11, R{x) shifts abruptly 
from one side of the domain wall to the other. Mermin|l| 
discusses some of the weird behavior one gets following 
this path. The point is, R{x) seems to be constrained: it 
doesn't change on its own, but follows the broken transla- 
tional order. Keeping it as an order parameter seems no 
more necessary than keeping p = j?/)! around in a super- 
conductor: only is massless, and p just wiggles around 
Pq in a boring way. 

Now, Ming and I have spent a huge amount of time 
trying to make these words into a mathematical theory. 
(We started with smectics, then studied superconduc- 
tors, then thought about some ideas of Toner and Nel- 
son, . . .) Ming has gone on to better things, and I'm 
still futzing with it. I can summarize where we are right 
now. Suppose we consider a rotationally distorted two- 
dimensional crystal (figure 12a). We can define a ro- 
tational order parameter by looking at the angle of the 
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nearest-neighbor bonds: 

i?(-)=f'°'' ""'V (15) 
\ — smO cos & J ^ ' 

The translational order parameter u is just as it always 
was: if x is the original position and p{x) is the corre- 
sponding position in the ideal lattice, 

u{x) — p{x) — X. (16) 

Now, the free energy can only depend on gradients of 
w, since it is translationally invariant. It also cannot 
change if we perform a uniform rotation: R RqR, 
p —^ RqP- From this, we can see that the free energy 
must be written in terms of gradients of R{x) and the 
particular combination p^sj 

2 

= Si-j - ^ Rktiduj/dxk + Skj). (17) 

k=l 

A reasonable free energy for a crystal then becomes 

Scrystal = (V^)^ + 2^ ^ (^^^^) ^ ^^^^^ 
ij 

This is just the normal elastic energy everybody uses, 
except for the third term multiplied by k. Normally, the 
strain matrix e is defined to be symmetric, so this term 
is then zero. 

Our free energy doesn't keep e automatically symmet- 
ric precisely because we have R{x) as an independent 
degree of freedom. The antisymmetric part measures the 
amount that R disagrees with the local gradients of H. It 
turns out that this antisymmetric part for the crystalline 
free energy is analogous to the current for the supercon- 
ductor, which has a Meissner effect. 

There are several things I haven't been able to do, 
though. First, I don't think ei2 — £21 is expelled quite like 
its analogue in the superconductor. I think we can show, 
though, that it is a boring variable like p was. Second, I 
haven't a clue on how to show that grains exist. To show 
that grains exist I have to show a constraint like = 0! 

We started this lecture by admiring the focal conic de- 
fects in smectic liquid crystals: beautiful ellipses and hy- 
perbolas which are due not to topology, but to geometri- 
cal consequences of a constraint. We saw how constraints 
can be enforced by the energy: "massive" modes decay 
exponentially. We saw explicitly how this occurs in su- 
perconductors — the magnetic field is constrained to zero 
because the photon and the Goldstone boson for the su- 
perconducting gauge symmetry combine into a massive 
particle. Finally, we discussed analogous effects in the 
everyday problem of grain boundaries in crystals, and 
realized that we don't really understand them in a deep 
sense. 



V. APPENDIX: THE SMECTIC ORDER 
PARAMETER 




FIG. 13: Equally Spaced Layers Imply curln = 0. Smec- 
tic layers, with a loop C enclosing an area A. The dot product 
h ■ dl gives the cosine of the angle of the curve C with respect 
to the layers, and a/ cos 6 is the length of curve C between two 
layers, so 1/a J^h-d£ gives the net number of layers crossed by 
the curve C. (A layer crossed first forward and then backward 
cancels, of course). Since in a closed loop the net number of 
layers crossed must be zero (assuming no dislocations), this 
must be zero. By Stokes' theorem, Jq'iT' ' d£ — curln • dA. 
This is true for any little area A, so curln = 0. 

Here we derive the consequences for layered systems of 
the constraint that the layers be equally spaced. Suppose 
that there are a stack of (bent) sheets, equally spaced 
from one to the next, with separation a. Suppose that 
the unit normal to these sheets at a position x is given by 
n. Consider traveling around a loop C, crossing various 
layers as we go around (figure Al). The number of layers 
we cross is given by the line integral 

(1/a) / n • d£ = net # crossed. (19) 
Jc 

If the layers exist throughout the region without any de- 
fects, then the net number crossed around any closed loop 
must be zero. Using Stokes' theorem, this integral over 
C is equal to an integral over the area A swept out by 
the curve: 

/ n-di^ / curln -dA. (20) 
Jc J A 

But for this to be true for all areas A, curln must be 
zero. 

Now, we already know that — 1. The derivative 
dfi^/dxa, of course, must be zero, so using the product 
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rule 

'^^fipdfip/dxa ^ 0. (21) 

Now, since we know curln = 0, we know from ^ that 

dfip/dxa = dfia/dxp. (22) 
Finally, combining these, we find 

'^^fipdfia/dxp — {h ■V)n = 0. (23) 

13 

This implies that n doesn't change when you move in the 
n direction. This means that fi will be the perpendicular 
to the next layer as well: that is, a straight line perpen- 
dicular to one layer will be perpendicular to every layer 
it crosses. 

These perpendicular lines are called generators. We 
qualitatively knew already that one layer determined its 
surroundings: now we have a simple geometrical rule de- 
scribing this nonlocal constraint. For your information, 
the defects occur where the generators cross (as shown 
in figure 3): this surface is called the evolute, or surface 
of centers, for the layer. 
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